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ABSTRACT 

Following  a  brief  discussion  on  the  background,  applicability, 
ar,J  limitations  of  th"  Vn Imogprov  Tort,  a  computer  program  of  the 
Kolmogorov  Test  for  normality  is  described.  The  program  is  applicable 
for  testing  the  null  hypothesis  that  a  random  sample  is  from  a  parent  normal 
population  whose  mean  and  variance  are  equal  to  those  of  the  sample 
distribution.  The  program  determines  the  minimum  and  maximum  sample 
values,  computes  sample  estimates  of  the  mean  and  variance  of  the 
hypothesized  normal  distribution,  and  cooputes  the  Kolmogorov  statistic. 

Five  optional  significance  levels  (0.20,  0.15,  0.10,  0.05,  end  0.01)  are 
available,  and  sample  size  limitations  are  4  £  n  -s  2500.  An  optional 
feature  provides  for  CRT  plot  output  applicable  for  both  hypothesis  testing 
and  interval  estimation. 

The  program  is  coded  in  FORTRAN  IV  for  the  IBM  7030  (STRETCH)  computer. 
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FOREWORD 

The  formulation  of  the  computer  program  for  the  Kolmogorov 
Test  was  performed  in  the  Mathematical  Statistics  Branch  of  the 
Operations  Research  Division.  The  coding  of  the  program  was  performed 
in  the  Operations  Sciences  Branch  of  the  Computer  Programming  Division. 
Dr.  H.  W.  Lilliefors  of  the  George  Washington  University  supplied  the 
critical  values  which  are  tabulated  within  the  program.  The  project 
was  supported  under  Foundational  Research  Project  No.  29Y,  "Computer 
Programs  for  Statistical  Analyses. " 

The  daLe  of  completion  was  15  December  1967. 
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I.  INTRODUCTION 


The  applied  statistician  is  often  confronted  with  the  problem  of 
investigating  the  distribution  of  sampled  data.  Because  of  the  desirable 
characteristics  of  the  normal  distribution  and  because  of  the  frequency 
of  occurrence  of  these  characteristics  in  physical  phenomena,  a  logical 
approach  to  the  problem  is  an  examination  of  the  sample  distribution 
to  determine  if  it  can  be  approximated  by  a  normal  distribution.  This 
examination  is  often  accomplished  by  hypothesis  testing,  l.e,,  testing 
the  null  hypothesis  that  the  randomly  observed  sample  date  la  from  a 
parent  normal  population. 

The  Chi-Square  Test  originally  proposed  by  Peerson  in  1900  eppeers 
to  be  the  most  popular  of  the  Goodness  of  Pit  Tests.  A  common  criticism, 
however,  of  the  Chi-Square  Test  ia  the  lots  of  power  because  of  the  requited 
grouping  of  the  sample  data,  e.g.,  see  Massey  (1951),  Birnbtum  (1952),  or 
Cochran  (1954)  .  In  addition,  to  employ  the  Chi-Square  Test  to  test  for 
normality,  a  minimum  of  four  intervals  (groups)  is  needed  when  the  population 
is  not  completely  specified  under  the  null  hypothesis  and  the  two  population 
parameters  (u  end  c*)  are  estimated  by  sample  statistics.  Consequently,  the 
Chi-Square  Test  is  not  applicable  to  very  small  samples. 

An  alternative  distribution-free  test  proposed  by  Kolmogorov  in  1933 
does  not  have  the  objectionable  requirements  mentioned  above  for  the  Chi- 
Square  Test.  That  is,  the  Kolmogorov  Test  requires  no  grouping  of  the  data 
and  the  test  applicable  to  very  small  samples,  Because  the  Kolmogorov 
Test  treats  individual  observations  rather  than  grouped  data,  sample 
information  may  be  better  utilized  by  the  Kolmogorov  Test  than  by  the 
Chi-Square  Test,  see  Blrnbaum  (1952).  Studies  by  Massey  (1951)  and 
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Kac ,  Kiefer,  and  Wolfowits  (1955)  indicate  that  the  asymptotic  power 
of  the  Kolmogorov  Test  is  greater  than  chat  of  the  Chi-Square  Test. 

A  recent  paper  fey  Slakter  (1967),  however,  shows  that  the  Kolmogorov  Test 
is  not  uniformly  more  powerful  than  the  Chi-Square  Test.  Further  studies 
are  necessary,  in  the  opinion  of  the  author  of  this  report,  to  Justify 
a  preference  of  either  test  on  the  basis  of  power.  On  the  other  hand, 
the  definite  advantage  of  the  Kolmogorov  Test  with  respect  to  small 
samples  is  sufficient  justification  for  the  applied  statistician  to  give 
serious  consiu. . -Li..;  to  the  Kolmogorov  Test. 

II.  THE  KOLMOGOROV  TEST 
A .  Background 

The  Kolmogorov  Test  for  normality,  like  the  Chi-Square  Test  for 
normality,  la  a  test  for  assessing  the  agreement  of  a  sample  distribution 
with  a  parent  normal  distribution.  While  the  Chl-Squarc  Test  is  concerned 
with  the  agreement  of  a  sample  distribution  with  a  theoretical  distribution, 
the  Kolmogorov  Test  is  concerned  with  the  agreement  of  a  sample  cumulative 
distribution  with  a  theoretical  cumulative  distribution. 

Consider  a  random  variable  X  with  a  continuous  cumulative 
distribution  function  (cdf)  ,  F(x)  ,  where 

x 

F(x)  =  P[X  s  x]  -  f  t ( x) dx  .  (1) 

•  00 

Let  Xj ,x? ,  ■  ■  ■ ,X)  , ■ •  ■  ,x  he  a  random  sample  of  n  observations  (or  measurements) 
of  X.  The  empirical  cdf  corresponding  to  the  random  sample  of  s tie  n  Is 
Sr (x)  ,  where 

Sr  (x)  =  number  of  x,  -  x 
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The  Kolmogorov  statistic,  D(n) ,  is  based  on  the  maximum  absolute  deviation 
between  F(x)  and  St (x)  .  That  is, 

D(n)  =  Max  j  F(x)  -  Sn(x),  .  (3) 

D(n)  then  is  the  test  statistic  to  test  the  nuLl  hypothesis  of  no 
difference  between  Sn(x)  and  F(x)  . 

Until  recently  the  Kolmogorov  Test  was  applicable  for  only  completely 
specified  F(x)  under  the  null  hypothesis.  That  is,  the  standard  tables, 
such  as  Massey  (1951),  Birnbaum  (1952),  and  Miller  (1956),  of  critical 
values  for  testing  D(n)  are  not  applicable  if  population  parameters  of 
the  theoretical  distribution  are  estimated  by  sample  statistics.  Lllliefors 
(1967)  provided  a  tabulation  of  critical  values  applicable  whan  the 
theoretical  distribution  is  normal  and  sample  statistics  are  used  to  estimate 
the  two  population  parameters.  The  critical  values,  which  were  deternined 
by  Monte  Car lo  calculations ,  are  for  five  significance  levels  (0.20,  0.15, 
0.10,  0.05,  and  0.01)  for  n  5s  4 .  Consequently,  the  Kolmogorov  Test  can 
now  be  used  to  examine  a  sample  distribution  to  determine  if  It  can  he 
approximated  by  a  normal  distribution  whose  mean  and  variance  are  equal 
to  those  of  the  sample  distribution. 
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B .  The  Kolmogorov  Teat  for  Normality 


To  distinguish  the  incompletely  specified  from  the  completely 
specified  theoretical  distribution  under  the  null  hypothesis,  a  circumflex 
is  used  denoting  the  use  of  sample  statistics  for  the  two  population  parameters 
of  the  normal  distribution.  Equation  (1)  then  becomes 

r(x)  =  P[X  s  x]  =  f  f(x)dx,  (A) 

•  0» 

where 

f<*>  -  (1//JW*)  a*p[-(«-K)a/2s2;),  (5) 

and  where  the  familiar  statistics,  X  and  s3,  are 

n  n 

X*  Ex,  /n  and  sa  ■  E(x, -^/(n-l)  .  (6) 

i«l  i«l 

a* 

Substituting  P(x)  for  P(x)  in  equation  (3),  the  Kolmogorov  statistic 
bee  otne  s 

D(n)  -  Max  |F(x)  -  SB(x)|  .  (7) 

The  null  hypothesis  that  a  random  sample  of  n  observations  is 
from  a  parent  normal  population  is  tested  by  comparing  D(n)  with  the 
critical  value,  Da(n),  of  TABUS  I.  If  D(n)  i  Da(n),  reject  the  null 
hypothesis  at  the  a-level  of  significance;  otherwise,  do  not  reject  the 
null  hypothesis , 


A 


The  Kolmogorov  Test  is  a  two-sided  test  in  that  departures  of 

A  * 

S„ (x)  in  either  direction  from  P(x)  increase  D(n) .  The  rejection 

A  A 

region  of  D(n)  ia  the  region  out aide  of  the  band  P(x)  +  Dc(n),  i  .e . , 


F(U)  =  F(x)  +  Da(n) 
F(L)  «  F(x)  -  Ds(n)  . 


(8) 


That  is,  the  null  hypothesis  is  rejected  at  the  o-level  of  significance 
if  S„(x)  passes  outside  this  band.  Because  Sn(x)  is  a  step-function, 
S„(x)  is  constant  for  x(4_j)  £  X  <  x(j),  where  x(j)  ia  the  jth  order 
statistic.  Therefore,  if  Sa(x)  passes  into  the  lower  rejection  region, 
D(n)  =  |F(x(j))  -  S,  (*( j-i))l  *Da(n),  see  Figure  1(a);  if  S„  (x)  passes 
into  the  upper  rejection  region,  6{n)  -  j F(x< jj)  -  Sn(x(J))|  2  D3(n),  see 
Figure  1(b) . 


Figure  . 
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In  addition  to  the  utility  of  Dfl{n>  for  hypothesis  testing,  Da(n) 
may  also  be  used  for  interval  estimation.  The  100(l-o,)7.  confidence 
limits  for  the  "true"  cumulative  distribution  function,  F(x) ,  are 

sa)  <  p(x)  <  s(u), 

where 

S(U)  =  S„(x)  +  Da(n) 

S(L)  =  S„(x)  -  Da(n> 

III.  COMPUTATIONAL  PROCEDURE 

Before  performing  the  computations  described  in  the  previous 
section,  the  n  sample  values  are  transformed  by  one  of  the  following 
thirteen  available  transformations. 


(9) 


(10) 


Transformation  Number 
1 
2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 


Transformation 
x  *-  x 
x  t  In  x 
x  «-  in(in(x)) 
x  *•  ln(A+x) 
x  *-  ln(B+ln(C+x) ) 
x  *-  /x 
x  *■  1/x 
x  ♦-  l/(D+x) 
x  «-  sin"'x 
x  ♦-  2  sin’Vx 
x  *-  x/E 
x  ♦-  sin  x 
x  *■  cos  x 
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The  eransf ormations  are  Identified  on  card  type  3,  and  the  constants 
A,B,C,D,  and  E  are  input  on  card  type  5  (see  Section  XV). 

After  ordering  the  n  transformed  sample  values,  the  subset 

*(1)  <  Hi)  <  <  H})  <  <  x(*>  :  k  s  0  (ll) 

of  sample  data  is  considered.  The  empirical  cdi,  Sn  (x) ,  is  evaluated 
for  each  of  the  k  unique  sample  values.  To  evaluate  F(x)  ,  the  k  unique 
sample  values  are  standardized  by  substituting  u  =  (x-7)/s  in  f(x) . 

The  standardized  pdf, 

$(u)  •  (L/m)  expC-Jpi8],  (12) 

aa 

is  then  used  to  evaluate  7(x)  which  becomes 

•  **  “ 

F(x)  =  *(u)  =  j  $( t)dt .  (13) 

mCO 

In  addition  to  X^jj,  x(*)»  x(*)"  x(l)i  n»  x>  sE  »  an(l  s*  system 

output  consists  of  u,  F(x) ,  S„(x),  and  |  F(x)  -  SB(x)|  for  each  of  the 
k  unique  order  statistics  along  with  the  identification  of  D(n)  .  CRT 
output  consists  of  one  or  more  of  the  four  following  plot  types. 

Plot  Type  A  -  SB(x),  F(x)  ,  F(U) ,  F(L) 

Plot  Type  B  -  SB(x),  F(x)  ,  S(U) ,  S(L) 

Plot  Type  C  -  SB(x),  S(U)  ,  S(L) 

Plot  Type  D  -  S„(x)  ,  F(x) 

The  desired  plot  type(s)  are  identified  on  card  type  3. 
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IV.  INPUT  PREPARATION 
A ,  Dtck  Setup 


The  input  deck  Is  listed  below  cy  card  type.  Multiple  jobs 
may  be  processed  by  stacking  card  types  I  through  6. 

CARD  TYK  I  -  JOB  HKNTXf  XCATION  CARD 
CARD  TYPE  2  -  VARIABLE  FORMAT  CARD 
CARD  TYPB  3  -  MAIN  CONTROL  CARD 
CARD  TYPE  4  -  ALPHA  IDENTIFICATION  CARD 
CARD  TYPE  3  -  TRANSFORMATION  CONSTANT  CARD 
CARD  TYPE  6  -  SAMPLE  DATA  CARD 

B .  Input  Deck  Dascription 

Entrlaa  for  varlablea  with  an  I  format  specification  must  be 
right-adjusted  in  the  specified  field. 


CARD  TYPE  1 

-  JOB  IDENTIFICATION  CARD 

Column 

Formet 

LHifikLf 

Description 

1-80 

1QA8 

JOI(I)  -  J0B(10) 

Job  Identification 

CARD  TYPH  2 

-  VARIABLE  FORMAT  CARD 

Column 

Formet 

Varleble 

Description 

1-80 

10AB 

FMT(l)  -  FMT(IO) 

Format  for  reading 
in  the  sample  data.  The 
format  specifications 
must  be  enclosed  in 
parentheses . 

CARD  TYPE  3 

-  MAIN  CONTROL  CARD 

Column 

Formet 

Verieble 

Descriotion 

1-2 

12 

NRUN 

Number  of  transformations 

to  be  performed 
<1  S  NRUN  <  13) 
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Column 

11 

12 


23 

25 

31 

32 

33 

34 

40-41 


Form* t  Variable 

n  XTRAN(l) 

r,  ITEAN(2) 


Dimer  lotion 

0  -  Do  not  perform 
transformation  number  1 
1  -  Do  perform 
transformation  number  1 

0  -  Do  not  perform 
transformation  number  2 
l  -  Do  perform 
transformation  number  2 


ITRAN  ( 13) 


MEDIUM 


12 


IA 


12 


a 


12 


IC 


12 


ID 


N0TIC 


0  -  Do  not  perform 
trana format ion  number  13 

1  -  Do  perform 

trans formation  number  13 

2  -  Sample  data  is 
input  on  cards 

1  -  Sample  data  la 
input  on  tape 

0  -  Do  not  perform 
plot  type  A 
1  -  Do  perform 
plot  type  A 

0  -  Do  not  perform 
plot  type  B 
1  -  Do  perform 
plot  type  B 

0  -  Do  not  perform 
plot  type  C 
1  -  Do  perform 
plot  type  C 

0  -  Do  not  perform 
plot  type  D 
1  -  Do  perform 
plot  type  D 

Number  of  desired  tick 
marks  on  the  abscissa 
axis  (if  left  blank 
N0TIC  is  set  i  qual  to  15) 
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Column 


Pnrnvu  > 


Variable 

IKA(l) 


46-47  X2 


IRA  (2) 


Humber  of  decimal 
digits  to  be  used  In 
labeling  abscissa  tick 
marks  for  transformation 
number  1 


Number  of  decimal 

digits  to  be  used  in 
labeling  abscissa  tick 
marks  for  transformation 
number  2 


68-69 

12 

ZltA(13) 

74-75 

12 

NBR 

78 

11 

NC 

Number  of  decimal 
digits  to  be  used  in 
labeling  abscissa  tick 
markB  for  transformation 
number  13 

Number  of  sample 
data  per  card 

0  -  Card  type  3  is  not 
input 

1  -  Card  type  j  is  input 


CARD  TYPE  4  -  AI.PHa  IDENTIFICATION  CARD 
To  obtain  plots,  at  least  one  alpha  must  be  indicated. 


SSEHt  variable 

11  IALPH4(1,I) 


IALPHA(2,1) 


Description 


1  -  Alpha  =  .20  Is 
desired  for 

transformation  number  1 
0  -  Alpha  *  .20  is 
not  desired  fur 
transforroation  number  1 


1  -  Alpha  -  .13  is 
desired  for 

transformation  number  1 
0  -  Alpha  =  .13  is 
not  desired  for 
transformation  number  1 
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Column 

3 


Format 

11 


Variable 

ULPmO.l) 


II 


lALPm(4,l) 


5  II  IALPHA<5,1) 


6  XI  IALPHA(1,2) 


Deacr ipt ion 

1  -  A lpha  ■  ,  10  is 
desired  for 

transformation  number  i 
0  •  Alpha  =  . 10  is 
not  desired  for 
transformation  number  I 

1  -  Alpha  ••  .03  is 
desired  for 

transformation  number  l 
0  -  Alpha  a  ,05  is 
not  desired  for 
trans format  ion  number  1 

1  -  Alpha  =  .01  is 
desired  for 

transformation  number  1 
0  -  Alpha  =  .01  is 
not  desired  for 
transformation  number  1 

1  -  Alpha  a  .20  is 
desired  for 

transformation  number  2 
0  -  Alpha  =  .20  is 
not  desired  for 
transformation  number  2 


65 


II 


IALPHA(5,13)  1  -Alpha  =  .01  is 

des  red  for 

tran:; formation  number  13 
0  -  Alpha  a  .01  is 
not  desired  for 
transformation  number  13 
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CARD  TVFE  5  -  TRANSFORMATION  CONSTANT  CARD 


This  card  type  is  omitted  if  NC  «■  0  on  card  type  3. 


Column 

Format 

Variable 

Descript  ion 

1-14 

F14.6 

A 

Constant  for 
transformation  number  4 

15-28 

F14.6 

B 

Constant  for 
transformation  number  5 

29-42 

F14.6 

C 

Constant  for 
transformation  number  5 

43-56 

F14.6 

D 

Constant  for 
transformation  number  8 

57-70 

F14.6 

I 

Constant  for 
transformation  number  11 

CARD  TYPE  6  - 

SAMPLE  DATA  CARD (S3 

The  sample 

size  (n)  must 

be  4  <  n  s  2500.  This 

card  type  is 

omitted  if  MEDIUM  =  1  on  card 

type  3. 

Column 

Format 

Variable 

Description 

1-2 

12 

ITEST 

Blank  -  This  is  not  the 
last  card  containing 
sample  data 

A)  -  Number  of  sample 
data  on  the  last 
card  or  record 

3-80 

Variable 

XX 

The  array  of  sample  data 

which  must  be  input 
according  to  the  variable 
format  on  card  type  2 
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c  •  Rsauest  Sheet* 


If  the  sample  date  (card  type  6)  ia  input  on  punched  card*,  the 


job  request  sheet  le  prepared  u  shewn  below. 
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2  1 

BD 

BB 

BB 

mnnn 

mfmm 

c  •  tr 
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T  A  1*  t  N  Uu  M  * 

Rffff 

h 

— 

MU  Motif.T  ax 

II 

MOODAHUUl  NUU|«  A 

1  1 

BUB 

SAICIAL  MAXOtlXfl 
/i*l  •! lie  In*  f . J 

hihold  r~iiouF  n *“  *«*»*•« 

*  0J  LJ4-0  LJilOl  fO»  ADO  I  P  t  DUAL 

joTciai  ihst»u”t  iomi  •»  tt¥trt») 
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If  the  temple  date  la  input  on  tape,  the  tape  number  is  written 
on  the  J*>b  request  sheet  instead  of  "Scratch".  And  tha  "RIEL,  NUL"  card 
(second  card)  of  the  IOD  deck  is  prepared  as  follows. 

Column 

1  Punch  tha  latter  "B". 

10-17  Punch  "UXL,  NUL"  (or  "UKL,  PUL"  if  the  tape  Is  file 

protected). 
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Punch  the  number  of  the  tape. 


The  CRT  request  sheet  Is  prepared  as  follows. 


TRAid  OUTPUT  MUM  OP  P »F(»  COPIES  BEQUEST 

2  iarv.  <U«At 


n  *OQH  AMMfcn 

ROOM 

BLDG. 

PHONE 

|  DA tt “ 

J  .  ORSUIAX 

A141 

1200 

73P.9-.  . 

L ttZ3Z£L_ 

P  ILM  IDENTIFICATION 


_ KS-A3 _ 

approximate  NUMar.n  or  frame* 

The  total  number  of  plots  requested  for  the  lobs  being  processed. 

number  of  paper  copies  per  frame 

The  nuaber  of  _gopie>  desired  of  each  plqt_._  _ 

FOR  OPERATORS  ONLY 

counter  reaoing  start  finish 


DATE  AND  TIME  PROBLEM  RAN 
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(Job  Identification  as  given  on  card  type  1) 
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KOLMOGOROV  TEST  FOR  NORMALITY 


Note:  f».\  =  F(.\  -  F(x/,\),  and  S 


B ,  CRT  Output 


The  plots  are  perforated  by  the  system  "GRF  Plot  Subroutine"  and 
are  plotted  on  the  CRT  printer  units .  The  plots  are  in  an  ordinary 
Cartesian  coordinate  system  with  origin  (x^  ,0)  .  The  abscissa  of 

the  system  is  the  x'-axis.  A  zero  (0)  is  plotted  at  the  point  (vt',0.5). 
The  four  types  of  plots  are  labeled  as  follows. 


Plot  Type  A  -  SB(x),  F(x)  ,  F(U),  F(L) 


OBSERVED  RELATIVE  CUMULATIVE  FREQUENCY, 

THEORETICAL  RELATIVE  CUMULATIVE  FREQUENCY,  AND 

CRITICAL  REGION  BOUNDS  FOR  ALPHA  -  a  LEVEL  OF  SIGNIFICANCE, 

D (ALPHA)  -  Da(n) 


Plot  Type  B  -  S„(x),  F(x)  ,  S(U),  S(L) 


OBSERVED  RELATIVE  CUMULATIVE  FREQUENCY, 

THEORETICAL  RELATIVE  CUMUIATIVE  FREQUENCY, 

100(1 -o)  PERCENT  CONFIDENCE  BANDS  FOR  THE  THEORETICAL 
CUMULATIVE  DISTRIBUTION  FUNCTION, 

D (ALPHA)  =  Da(n) 


Plot  Type  C  -  S,(x),  S(U),  S(L) 


OBSERVED  RELATIVE  CUMULATIVE  FREQUENCY  AND  THE 
100(1 -a)  PERCENT  CONFIDENCE  BANDS  FOR  THE  THEORETICAL 
CUMULATIVE  DISTRIBUTION  FUNCTION 


Plot  Type  D  -  S„(x),  F(x) 

OBSERVED  AND  THEORETICAL  RELATIVE  CUMULATIVE  FREQUENCY 
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The  running  time  1b  printed  *t  the  end  of  each  job  processed.  The 


following  table  of  observed  times  may  be  used  as  a  guide  for  estimating 


program  running  times. 

Sample  Size  No.  of  Runs  Total 

No.  of  Plots 

Time  in  Seconds 

25 

1 

4 

12  .4 

140 

1 

4 

16.1 

10 

1 

6 

19.0 

30 

1 

6 

20.0 

20 

2 

8 

25.4 

20 

4 

16 

45.8 

20 

5 

20 

56.4 

1000 

1 

2 

54.6 

VI.  EXAMPLE 

PROBLEM 

A.  Problem  Description 

The  example  problem  consists  of  20  observations  (x(  -values)  randomly 
sampled  from  a  continuous  distribution.  Two  runs  were  performed,  one  on  the 
original  observations  and  one  on  the  logarithmically  transformed  values. 
Because  the  minimum  original  sample  value  was  -750,  transformation  number  4 
was  used  with  the  constant  A  =  751,  All  four  available  plot  types  were 
performed.  The  0.05-level  of  significance  was  chosen  for  testing  the  null 
hypothesis  of  normality. 

The  input  deck  setup  is  illustrated  on  the  following  Data  Card 
Layout  Sheet. 
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Program  Output 
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£14MPI£  PBOH.C* 

RUN  Hy*»e£«  a  WITH  20  OBSCBVAl  1  OHS 

0»S£B»t9  ANO  TmMBETICAL  BELATIVE  CUMULATIVE  fRMutMCT 
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T  1MK  IN  SECONDS  = 


ICO 


RUM  MIMIC*  2  WITH  20  I0W 

OIS.CRVCD  *8.  At  I  VC  com*.  AT  I  *C  FRCOUCNCV  AND  T  HC 

«  *C*CCMT  COMCiCCMCC  BANDS  FO*  THC  THCORCTICM.  CI*H*.Atl*C  BIST* 


ci  ample  Pft0M.cn 

■UN  MWICI  2  WITH  20  OtSEBVAl IONS 

0»SE»v£0  AND  IMtOOCTICAL  BCLAllVC  CUNU.Al|«C  f«Ml<NC> 


D .  Discussion  of  Teat  Resales 

To  test  the  null  hypothesis  “The  20  randomly  sampled  observations 
are  from  a  parent  normal  population"  at  the  0.05-level  of  significance, 

D(n)  b  0.2568  (from  the  first  run)  is  compared  with  D0  ...,(20)  -  0.190.  The 
null  hypothesis  is,  therefore,  rejected  at  the  0.05*level  of  significance. 
Note  that  the  maximum  absolute  deviation  appears  at  X(10)  =  -699,  and  at 

this  point  Sa(x)  has  passed  above  F(U)  into  the  upper  rejection  region 
(see  pages  24  and  25) .  To  test  the  null  hypothesis  of  normality  of  the 
transformed  random  variable  at  the  0. 05-level  of  significance,  D(n)  =  0.1615 
(from  the  second  run)  is  compared  with  Do>ob(20)  =  0.190  (see  page  31); 
consequently,  the  null  hypothesis  is  not  rejected.  Unlike  run  number  !, 
plot  type  A  for  run  number  2  on  page  32  shows  S„(x)  remaining  within  the 
"acceptance"  band.  Therefore,  we  conclude  at  the  0.05-level  of  significance 


that  the  original  random  variable  is  approximately  log  normally  distributed. 
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VII.  PROGRAM  LISTING 


T YPE. COMP l LCO. FORTRAN, PM 
SUBTvPE.FlOO 

10  1 00* TAPE*  « ,  .EVEN*  .SAVE 

REEL .NULXXXX 
2  I OO. •READER 

ENO 

SUBTYPE • FORTRAN «  LMAP  *  PUNCH 
KOLMOGOROV  TEST  FOR  NORMALITY 
DIMENSION  TABLE (28*  5) 

DIMENSION  JOBdO),  FMT  (1 0 1  «  IFILEOB),  1TRANU3).  I  ALPHA  (S,|3>* 

1  X  (2500  )  <  U  <  2500  >  •  ALPHA  <51*  1 FREQ (2500),  TEN(20>*  TC2S01)*  IRAI13 

2  ).  XX ( 2500 ) «  S ( 2500 ) «  SNK2500),  SN(2500>,  F(2900>»  FN(SOO>. 

3  FL(500).  FU(500 ) *  SL(2S00),  SU<*500) 

COMMON  /PRIME/  XP (2*  1290).  XXP(2*  12501 
EQUIVALENCE  (IFILE(l),  NRUN)«  (IFILEI2I*  I  TRAN  1 1)1* 

1  (  IFILE (3)*  I  TRAN ( 2  )  ) ,  (IFILEI4)*  1TRAN<3>),  (IFILE(S)*  |TRAN(4))» 

2  I  IF 1LE (6  t •  I  TRAN (9 )  )  ,  <IFILE(7>.  I  TRAN (6  I ) *  (1FILE<8>,  ITRANI7>>« 

3  (  IF  I  LE  (9  )  •  I  TRAN  (8)1*  (1FILEII0I*  I  TRAN  (9  I  >  «  (IFILEUI).  ITRANUO 

4  ))*  I  IFILE  (  )2  )  *  ITRAN(U))*  (  IFILE  (13)*  I  TRAN  ( 12  ) )  ,  (IFILE(|4). 

5  I  TRAN  (  13  )  )  *  (  (FILE  (  1  S  I  «  MEDIUM)*  (1FILEU6)*  1A)«  (IFILEU7)*  IB) 

6  •  (IFILE(IB),  IC),  (IFILEU9).  ID)*  IIFILE(20)*  N0T1C) 

EQUIVALENCE  (IFILE(2t)«  IRA  ( 1  )  )  •  (IFILE(22>.  IRA(2M*  <IFILE<23)« 

1  IRA (31)*  ( I F I LE ( 24 ) •  IRA(4))«  (IF|LE(29)«  |RA(5))«  MF1LE(26>, 

2  IRA (6  > ) •  ( IFILE (27) *  1RA<7)>,  (IFILE(2B>«  IRA(8)>*  <»F|LE(29). 

3  IRA  (9)  )  «  (IFILEOO).  IRA  (10))*  I  IFILE  (31  )  *  IRA(ll))*  <!FILE<32). 

4 IRA( 12 ) ) •  ( IFILE (33) •  tRA(13))«  (IFILE(34)«  NSR ) *  (IF1LE09).  NC) 

equivalence  (xp  ( i«i)«  x  ( i ) )  •  rxxPd.n*  xxu)) 


DATA  I  ALPHA ( I ) .  I  - 

I,  S) 

(.20, 

•  19, 

•  10, 

.05*  .01) 

DATA  (TABLE (1*1)*  I 

•  1 , 

28) 

11.300.  .285,  .265. 

•247, 

•  233. 

*223, 

•  219* 

•206, 

•  199. 

•  190, 

•  183* 

2  *177.  .173*  .169* 

•  166, 

•  163. 

•  160. 

.160. 

•  160, 

•  160. 

•  160, 

•  142* 

3  *142.  *142*  .142, 

.142, 

.  131  , 

•  736) 

DATA  (TABLE  (1*2). 

1*1. 

28) 

1I.319,  .299,  .277* 

.258, 

•244* 

•  233, 

•224* 

•  217, 

•  212. 

•202, 

•  194* 

2  *187,  .182,  .177* 

.173, 

•  169, 

•  166, 

•  166, 

•  166, 

•  166. 

•  166, 

*147* 

3  .147,  .147,  .147, 

.147. 

•  136, 

•  768) 

DATA  (TABLE  (1.3). 

I  ■  1 , 

28) 

1 ( *  352 ,  .3)5,  .294, 

*276, 

•  261. 

•  249, 

•239, 

•  230  * 

•  223. 

•  214, 

.207, 

2  .201 .  *195*  .189* 

.104, 

•  179. 

•  174, 

•  174. 

•  174, 

•  174. 

.174, 

•  158, 

3  *158.  *158.  .158* 

.198, 

•  144, 

•  805) 

DATA  (TABLE  (1,4), 

I  ■  1  , 

28: 

1 (.381 ,  .337*  .319, 

.300, 

.289* 

.271  , 

•258, 

•  249, 

•  242, 

•234, 

.227, 

2  .220,  .213,  .206, 

.200, 

•  195, 

•  190, 

•  190, 

*190. 

•  190, 

•  190, 

•  173, 

3  .173.  .173,  .173, 

.173, 

•  161  , 

*886) 

DATA  (TABLE  (1.5). 

1  •  1  , 

28) 

1  (.*'17,  .405,  .364, 

.348. 

•  33 1 , 

.3)  1  • 

•294, 

.284, 

•  275. 

•266, 

*261 , 

2  *257,  .250*  *245, 

.239, 

•  239, 

*231  , 

.231 . 

•  231  , 

•  231 . 

•  231, 

•  200, 

3  .200,  .200,  .200, 

.200, 

.187, 

1.031  ) 

CALL  CRT  ID  (4HKSA3) 

CALL  GRS  t  (.06652, 

.06652 

,  1  • , 

.06652) 

CALL  SET  IT 

CALL  MOVE  (0.,  XP, 

5000  ) 

READ  INPUT 

too  READ  no,  JOB,  FMT 

110  FORMAT  (10A8/  10A6) 

READ  120,  IFILE 

120  FORMAT  (12,  BX,  1311,  IX, 

11,  5X«  4)1 

.  5X, 

12,  2X«  1312,  4X« 

12. 

1 2X «  It) 
READ  130, 


I  ALPHA 
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130  FORMAT  <6511 > 

IF  <NC )  160*  160*  160 

160  READ  150.  A,  6*  C*  O.  E 
1 50  FORMAT  ( 5F 16*6) 

|60  K  ■  I 

KK  *  NCR 

|70  READ  (MEDIUM*  FMT )  ITEST*  < XX <  I  1  *  1  ■  K*KK) 

180  IF  (1TESTS  200*  190*  200 
!90  K  «  K  +  NBB 
KK  *  KK  +  NBR 
00  TO  170 

200  N  m  K  +  ITEST  -1 
210  1  RUN  *  1 

C  TRANSFORMATION  1 

IF  (ITRANdl)  220*  260*  220 
22 0  DO  230  I  -  1  *  N 

230  X( I  I  «  XX < I ) 

IT  «  1 

ASSIGN  260  TO  ILINE 
GO  TO  740 

C  TRANSFORMATION  2 

240  IF  < I TRAN (2)1  250*  29 0*  2S0 
250  IT  -  2 

ASSIGN  290  TO  ILINE 
DO  280  I  -  1 «  N 
IFIXXIIII  260*  260*  280 
260  PRINT  270*  IT 

270  FORMAT  (24H  ERROR  IN  TRANSFORMATION,  3X.  12*  3X.21HLN  OF  NEGATIV 

IE  NUMBER) 

GO  TO  1330 

280  X< I  >  ■  ALOG (XX  I  I ) ) 

GO  TO  740 

C  TRANSFORMATION  3 

290  IF  ( 1  TRAN 1 3 ) )  300,  330*  300 
300  IT  ■  3 

ASSIGN  330  TO  ILINE 
DO  320  I  ■  1  *  N 
IF  <XX< II)  310*  320.  320 
310  PRINT  270*  IT 
GO  TO  1330 

320  XII)  -  ALOG (  ALOG < XX ( I ) ) ) 

GO  TO  740 

C  TRANSFORMATION  6 

330  IF  ( I  TRAN ( 4 ) )  340*  370*  340 
340  IT  ■  4 

ASSIGN  370  TO  ILINE 
DO  360  I  ■  1  *  N 
AR6  «  A  +  XX ( l ) 

IF  (ARG)  350,  360*  360 
350  PRINT  270*  IT 
GO  TO  1330 

360  XII)  >  ALOG  (ARG) 

GO  TO  760 

C  TRANSFORMATION  5 

370  IF  ( 1  TRAN ( 5 ) )  380*  620*  380 
380  IT  «  5 

ASSIGN  620  TO  ILINE 
DO  40C  I  ■  1  ,  N 
ARG  «  C  ♦  XX ( I ) 
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390 


IF  l ARG)  410*  390*  390 
ARG2  ■  B  +  AL0G4ARG) 

IF  <ARG2)  4)0*  400*  400 
400  XI  I  I  -  ALOG  ( ARG2  > 

GO  TO  740 

410  PRINT  270*  IT 
GO  TO  1330 

C  TRANSFORMATION  6 

420  IF  t I  TRAN (6) I  430*  470*  430 
430  IT  •  6 

ASSIGN  470  TO  ILINE 
00  460  I  ■  I «  N 

IF  (XXtlU  440  *  460  *  460 
440  PRINT  450*  IT 

450  FORMAT  <24H  ERROR  IN  TRANSFORMATION*  3X«  12*  3X .  32HSQUARE  ROOT 
I  OF  A  NEGATIVE  NUMBER) 

GO  TO  1330 

460  XI  I )  >  SORT  I XX 1 1 ) ) 

GO  TO  740 

C  TRANSFORMATION  7 

470  IF  ( I  TRAN 17))  480*  520*  400 
400  IT  ■  7 

ASSIGN  520  TO  ILfNE 
00  S10  I  •  I*  N 

IF  ( XX ( I ) )  510*  490*  510 
490  PRINT  500*  IT 

500  FORMAT  (24H  ERROR  IN  TRANSFORMATION*  3X*  12*  3X«  1IHZERO  DIVIDE) 
GO  TO  1330 

510  XI  I )  ■  1*0/  XXtl ) 

GO  TO  740 

C  TRANSFORMATION  8 

520  IF  ( I  TRAN 1 8 ) )  530.  560.  530 
530  IT  ■  B 

ASSIGN  560  TO  ILINE 
DO  550  I  ■  I,  N 
ARG  ■  D  +  XXI 1 > 

IF  (ARG)  550*  540*  550 
540  PRINT  500*  IT 
GO  TO  1330 

550  XII )  >  1*0/  ARG 
GO  TO  740 

0  TRANSFORMATION  9 

560  IF  ( 1  TRAN 1 9 ) )  570.  610*  570 
570  IT  ■  9 

ASSIGN  610  TO  ILINE 
00  600  I  «  I «  N 

ARG  «  1 •  -  XX  < I  )  **  2 
IF  (ARG)  580*  580*  600 
580  PRINT  590.  IT 

590  FORMAT  (  24H  ERROR  IN  TRANSFORMATION.  3X *  12  *  3X.  47HZER0  DIVIDE 

I  OR  SQUARE  ROOT  OF  A  NEGATIVE  NUMBER) 

GO  0  1330 

600  X  < I  )  ■  ATAN  (XXII)/  SORT  (ARG)) 

GO  TO  740 

C  TRANSFORMATION  10 

610  IF  ( I  TRAN  1 1  0 ) >  620*  650*  620 
620  I T  ■  1 0 

ASSIGN  650  TO  ILINE 
DO  640  1  ■  1 •  N 
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ARC  •  I •  -  XX  (  I  ) 

IX  (ARC  ,GT*  0«  *  AND*  XXCi  *G£*  0.)  CO  TO  6*0 
630  PRINT  590*  IT 
CO  TO  1330 

6*0  X(1)  ■  2*  #  ATAN  (SORT  (XX (I ) I/SORT (ARC > ) 

CO  TO  7*0 

C  TRANSFORMATION  II 

6S0  IF  ( I  TRAN (It))  660*  690*  660 

660  1 T  »  1 1 

ASSIGN  690  TO  ILINE 
00  660  I  ■  I  *  N 
IF  (E>  660*  670.  650 
670  PRINT  BOO*  IT 
CO  TO  1330 
660  X « t  »  *  XX C I  )/E 
CO  TO  7*0 

C  TRANSFORMATION  12 

690  IF  ( I TRAN( 12)1  700*  720*  700 

700  DO  710  I  ■  I .  N 
710  X  <  I  >  ■  S I N  ( XX  (  I  )  ) 

IT  ■  12 

ASSIGN  720  TO  ILINE 
GO  TO  7*0 

C  TRANSFORMATION  13 

720  00  730  I  ■  1 ,  N 

730  X< I  I  ■  C0S(XX( I ) 1 
IT  »  13 

C  PRINT  IWHJT 

7*0  PRINT  750*  JOB*  FMT 

750  FORMAT  (1M1*  10A8/  23H  FORMAT  FOR  SAMPLE  DATA*  5X.  10*6) 

PRINT  760*  NRUN*  IRUN*  (IFILE(I)*  I*  2*  19) 

760  FORMAT  (15H  NO*  OF  RUNS  ■  «  12*  5X.  7HIRUN  a  *  12/ 

121H  TRANSFORMATION  |  »  «  II*  5X.  20MTRANSFORMAT I  ON  2  «  *  II*  5X, 

2  20HTRANSFORMATION  3  ■  ,  II*  5X»  20HTRANSFORMAT I  ON  4  *  ,  11/ 

3?1H  TRANSFORMATION  5  ■  «  II*  BX«  20HTRANSF0RMAT I  ON  6  ■  .  II*  5X* 

*  20MTRANSF0RMAT I ON  7  ■  •  II*  SX*  20HTRANSF0RMAT I  ON  8  *  .  11/ 

521 H  TRANSFORMATION  9  »  «  II.  5X*  20HTRANSF0RMAT I ON  10  *  .  II*  5X« 

6  20HTR ANSFORMAT I  ON  11  a  ,  M«  5X*  20HTRANSFORMAT |ON  12  "  .  11/ 

721 H  TRANSFORMATION  13  a  ,  It*  5X.  9HME0IUM  »  ,  11/ 

61BH  PLOT  TYPE  A  a  ,  1 1,1  lx.  1 AHPLOT  TYPE  B  »  *  I  1 « 1 1 X •  14HPL0T  TYP 

9E  C  a  ,  I  |,1  IX*  1  AHPLOT  TYPE  Da  ,  M) 

PRINT  770,  (IFILE(I),  I  a  20*  35) »  I  TEST,  I  ALPHA 
770  FORMAT  <  9H  NOTIC  a  ,  12.  1 6X *  6HIRA  •  ,  13(  12.  IX)/  7h  NBR  «  , 
1I2.16X.  5HNC  a  *  I),  20X *  BHI TEST  «  *  I2/24H  ALPHA  IDENTIFICATION 
2a  ,  13(511*  IX)) 

IF  <NC)  760*  800*  780 
760  PRINT  790,  A,  8.  C,  0,  E 

790  FORMAT  (29H  TRANSFORMATION  CONSTANT  CARD*  3X*  AHA  a  ,  FI  4. 6.  3X. 

1  4H0  a  .  F14.6,  3X*  4HC  a  ,  F14.6/32X.  4H0  «  •  F14.6,  3X ,  4HE  =  . 

2  F14.6) 

600  NM  a  N/2  +  MOD  (N*2) 

NN  a  0 

eiO  NX  a  NN  +  1 

NN  a  Ml  NO  (50.  NM )  +NN 
C  PRINT  ORIGINAL  AND  TRANSFORMED  DAT* 

PRINT  620, ( (XXP(J*I I ,  XP(J,I),  J  a  1,2),  I  «NX , NN ) 

820  FORMAT  ( 1 H 1 ,  9X ,  1 3H0R1G I NAL  DATA,  15X,  I 6H TRANSFORMED  DATA.  )6x. 
1I3HORICINAI  DATA,  15X,  1 6HTRANSF0RMED  DATA/  (SX,  F20.10.  10X,  F2 o. 

210,  I0X,  F20.10,  IOX,  F20.10I) 
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m  •  m  -  so 

if  (NM)  339*  830,  810 
L  SORT  DATA 

630  K  ■  N  -  1 

DO  850  I  «  1 ,  K 
ll  •  l  +  l 
DO  650  J  m  ll,  N 
1 F  (X(ll  -  XU>)  850,  850,  8*0 
640  TEMP  ■  X(J) 

X( J)  *  X< I ) 

X( ! »  *  TEMP 
850  CONTINUE 

C  ESTABLISH  PARAMETERS  FOR  PLOTS 

XM I  *  XU) 

XMA  «  X  (N ) 

RANGE  •  XMA  -  XM I 

IF  <  NOT  I C  ,EQ.  0>  NOTIC  •  15 

OX  «  RANGE  /  FLOAT (NOTIC) 

TEN  (1)  »  tO.E-IO 

IF  <DX  *OT«  TEN ( I  )  )  GO  TO  860 

OX  «  TENO  ) 

GO  TO  900 

860  IF  (DX  ,LT.  10,E10>  GO  TO  870 
OX  ■  10.E10 
GO  TO  900 

870  DO  880  M  ■  2,  20 

TEN  (M)  >  10.**  (M  -  m 

IF  (OX  rGT •  TEN(M-I)  .ANO.OX  *LE.  TEN (M )  )  GO  TO  890 
880  CONTINUE 

890  NOX  ■  DX/  TEN (M—l )  +  ,9 

DX  «  FLOAT  (NOX)  *  TEN(M-l) 

900  DELTA  *  RANGE/  500. 

FLN  ■  FLOAT (N) 

SUM  -  0. 

SUMSO  «  0* 

C  COMPUTE  AND  PRINT  MEAN  ADO  VARIANCE 

00  910  I  «  I ,  N 
SUM  ■  SUM  +  X ( I ) 

9 1 0  SUMSO  ■  SUMSO  +  X  (  I  )  *  XU) 

XBAR  »  SUM/  FLN 

SIGMAS  ■  (SUMSO  “  FLN  *  XBAR  *#  2)/  (FLN  -  1.) 

SIGMA  •  SORT  (SIGMA2) 

PRINT  920,  XM! ,  XMA,  RANGE,  N,  XBAR,  SIGMA2,  SIGMA 
920  FORMAT  (1H0,  BX«  1QHMINIMUM  •  ,  E15.8,  20X,  10HMAXIMUM  *  ,  El5,8« 
1I0X,  8HRANGE  «  ,  EJ5.8/  5X,  14HSAMPLE  SIZE  *  ,  15,  33X,  7HMEAN  » 
2,  E 1 5 , 8/  8X»  ItHVARIANCE  »  ,  E15.8,  9X,  21 HSTANDARD  DEVI  ATI  IN  a  , 
3E1S.8) 

C  COMPUTE  EMPIRICAL  CDF 

DO  »  0, 

J  ■  I 

I FREQ  <JS  ■  0 
DO  990  I  «  1  ,  N 
IFREO(J)  ■  IFREO(J)  ♦  I 
IF  (X(I)  -  XU+1  ))  930,  990,  930 
930  XU)  ■  X<  1  ) 

UU)  ■  (X(J  -  XBAR  >/  SIGMA 

F(J)  «  FREO  (Uu)  f 

SU)  *  FLOAT  (I)/  FLOAT  (N) 

C  DETERMINE  TEST  S TAT  IS  (C  D 
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SN(  J )  ■  ASS  (F( J)  -  hi  J) ) 


IF  (J  -  1  )  950,  990.  950 
940  SN1  (I )  «  ABS  (F{ J) 

60  TO  960 

950  SN1(J)  »  ABS  (F(J)  S(J  “III 
960  ZMAX  «  AMAX1  (SN(j;  SNltJ)) 
IF  (ZMAX  -  00)  9B0  >80,  970 

970  JJ  *  J 

00  •  ZMAX 
980  J  ■  J  +  1 

IF9E0  (J)  •  0 
990  CONTINUE 
C  PRINT  STATISTICS 


J  •  J  -  l 
II  «  J 
IK  *  0 

1000  IL  «  IK  +  | 

IK  «  M l NO  (50.  II)  +  IK 
PRINT  1010 

1010  FORMAT  (1MI,  59X.  29HKOLMQGOROV  TEST  FOR  NORMAL  I TY//3X .  IIHOR 

IOER  INDEX.  6X.  1 1 HTRANSF ORMEO . 1 8X •  1 2HSTANDAR0I ZED .  IX. 

2  I5HTHE0&ETICAL  CDF,  2X.  12H0BSERVED  CDF*  3X.  14HA8S0LUT 

3E  VALUE.  3X«  14H ABSOLUTE  VALUE/  7X,  3H(J»,  11X,  9HDATA  X(J>. 

4  7X,  9HFRE0UENCY.  4X. 

5I2HVARIATE  U(J)«  5X.  4MF(J),  12X.  4HS(J).  9X.  II HF ( J )-S ( J- I > ,  7X, 
69HF(J)-S(J)/I 

PRINT  1020. <  < I »  X<1).  I FREQ  (I),  U(I).  F(I),  S(I).  SN1  ( I ) ,  SN( I ) 

I).  I  -IL.IK) 

1020  FORMAT  (6X,  14,  6X.  F14.6.  I0X.  14,  4X.  Fl4.ll,  IX.  F14.12, 

I  2X«  F14.12,  IX,  F14.13,  2X,  F14.13) 

II  -  II  -  50 

IF  (II)  1030,  1030,  1000 
1030  PRINT  1040,  JJ,  DO 

1040  FORMAT  t  3H0D ( ,  14,  4H )  «  ,  F7.4) 

T ( l )  «  XMI 

C  COMPUTE  THEORETICAL  CDF 

DO  1050  L  ■  1.  500 
UV  *  (T (L )  -  XBAR ) /S I GMA 
FN(L>  ■  FREQ (UV ) 

IF  (T (L  >  *GT •  XMA)  60  TO  1060 
T  (L+l  )  •  TIL)  OELTA 
1050  CONTINUE 
1060  DO  1320  II  «  1.9 
C  OBTAIN  VALUE  FOR  D  SUB  ALPHA  OF  N 

IF  (I ALPHA  (It,  IT))  1070,  1320,  1070 
1070  IF  (IA  »EQ.  0  .AND.  IB  .EQ.  0  .AND.  IC  ,EQ.  0)  60  TO  1150 
IF  (N  -  30)  1080,  1080,  1090 


1080  NN  «  N  -  3 


OALPHA  >  (TABLE  (NN,  II)) 

GO  TO  1100 

1090  OALPHA  ■  TABLE  (28,  II)/  SORT (N ) 

1100  IF  (IB  .EQ.  0  .AND,  IC  .EQ.  0)  GO  TO  1120 
C  COMPUTE  CONFIDENCE  LIMITS 


DO  1110 
SU ( I  )  ■ 
1110  SL ( t  )  a 
1120  IF  (1A) 
C  COMPUTE 

1130  DO  1140 


I  ■  1,  J 
S ( I )  +  OALPHA 
S(  I )  -  OALPHA 
1130,  1150,  1130 
REJECTION  REGION 
I  «  1 ,  500 
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o  u 


FU(I t  «FN(|>  +  DALPHA 


1140  FLU)  ■FNIII  -  OALPHA 
**  *#*P  LOTS 
HEADINGS 

1 | SO  OC  1310  XX  «  1.4 

IP  (1PILE  (KK  +  1511  1160.  1310*  1160 
1160  PAINT  1170.  JOB  t  IRON,  N 

1170  FORMAT  (2HP2.  |5X,  I0A8/  1HP,  |SX.  UHRON  NUMBER  .  12.  6M  WITH  « 
!!5.  I 3H  OBSERVATIONS! 

GO  TO  < 1  ISO.  1200.  1220.  1240>.  KK 
C  HEADINGS  FOR  PLOT  TYPE  A 

1180  PRINT  1190.  ALPHA! IIJ,  DALPHA 

1190  FORMAT  (IHS.15X.  39H0SSCRVED  RELATIVE  CUMULATIVE  FREQUENCY  »/  1HP. 
USX,  40HTHE0RET I  CAL  RELATIVE  CUMULATIVE  FREQUENCY.  AND/  |HP, 

2 1 5X,  34HCRI T I CAL  REGION  BOUNDS  FOR  ALPHA  ■  »  PS. 2.  3X. 

^  22HLEVEL  OF  S IGN IF ICANCE ./I  HP,  15X.  11H0IALPMA1  •  ,  FS.3) 

GO  TO  1260 

C  HEADINGS  FOR  PLOT  TYPE  B 

1200  ICB  «  100*  •  (1.  -  ALPHA! Ill)  ♦  .3 
PRINT  1?10<  ICB,  DALPHA 

1210  FORMAT  (1HP.1SX.  39H0BSERVE0  RELATIVE  CUMULATIVE  FREQUENCY,/  |HS, 
USX,  42HTHE0RET1CAL  RELATIVE  CUMULATIVE  FREQUENCY./  1HP. 

2ISX.  12. 1X.TBHPERCENT  CONFIDENCE  BANDS  FOP  THE  THEORETICAL  CUMULAT 

aivE  distribution  function,/  ihp.isx.  iihdialpha)  «  .  fs.3> 

GO  TO  1260 

c  headings  for  plot  type  c 

1220  ICB  -  100.  *11.-  ALPHA  !l  in  *  .8 
PRINT  1230.  ICB 

1230  FORMAT  (IHP.ISX.  46MO0SERVED  RELATIVE  CUMULATIVE  FREOUENCY  AND  THE 
1/  IHP.1SX.12,  IX. 77HPEPCENT  CONFIDENCE  BANDS  FOR  THE  THEORETICAL  C 
2UMULATIVE  DISTRIBUTION  FUNCTION) 

GO  TO  1260 

C  HEADINGS  FOR  PLOT  TYPE  O 

1240  PRINT  1250 

1250  FORMAT  (  1HX.1SX,  54H0BSERVE0  AND  THEORETICAL  RELATIVE  CUMULATIVE 


1 FREOUENCY ) 

C  PLOT  MEAN  AND  S(X 

1260  CALL  GRF  !XMl,  0. 
I  .5,  1,  1,  2H 12, 
CALL  GRF  (XMJ,  0. 
1  S ( 1 ) ,  1 «  1 «  2H7S 
GO  TO  (1270.  12B0 
C  ADDITIONAL  CURVES 

1270  CALL  GRF  (XM! ,  0. 
1FNU),  1,  1,  2H75 
CALL  GRF  (XM1,  0. 
1FLU).  1.  1.  2H75 
CALL  GRF  ! XM 1 .  0. 


XMA,  ]•«  XM I ,  0*. OX, .OS .IRA (IT). 2.0.  1,  XBAR, 

1.0) 

XMA.  !.,  XMJ,  0 • «  DX. .05. IRA! IT) ,2. 1 ,  J,  X(l), 

0,  0) 

1290,  1300),  KK 
FOR  PLOT  TYPE  A 

XMA.  |«,  XM 1 ,  0*.  DX, ,05. IRA! IT) ,2, 1 . 1000. T  (1  )  ♦ 

-1  .  0) 

XMA,  I.,  XM I ,  o.«  OX..0S. IRA! IT) ,2. 1 « 1000.T  II  ). 

-1  .  0) 

XMA.  1,.  XM1,  0.,  DX, ,05. IRA! IT). 2.1,1000. Til )• 


1FUII).  1.  1.  2H75 
GO  TO  1310 

C  ADDITIONAL  CURVES 

1260  CALL  GRF  (XM] ,  0. 
1FN! 1 ) ,  1 .  1 .  2H75 
CALL  GRF  ( XMJ ,  0. 
1 SL (11.  1 .  1 .  2H75 

CALL  GRF  (XMI,  0. 
1 SU( 1 ) .  I .  1 ,  2H75 
GO  TO  1310 

C  ADDITIONAL  CURVES 


*  i  »  u  i 


FOR  PLOT  TYPE  B 


XMA, 

-1  «  0) 

1  «  « 

0) 

XMA. 

0,  0) 


; xm 1 1 


XMI  , 


XMI 


FOR  PLOT  TYPE  C 


DX.. 05. IRA (IT) .2,1 • 1000.T ( 1 > . 
DX, .05. IRA! IT) ,2, 1 ,  J,  X(l). 
DX,  .05.  IRAUT)  .2.  1  ,  J,  XU). 
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XM1  .  0.*  DX  •  •  OS  *  I  RA  (IT)«2tlt  J.  XU)* 


1290  CALL  GRF  (XM(*  0.,  XMA.  1 
1 SL  1 1 ) «  I*  1*  2H73 «  0,  0) 

CALL  GRF  <XMt  *  0.  .  XMA  «  !•*  XM]«  0»«  OX  ♦  .05 .  IRA  <  I  T  >  ,  2.  1  «  J.  XU)* 
ISUC1  >  ♦  I  .  1  »  2H75 1  0*  0) 

GO  TO  1310 

C  ADDITIONAL  CURVES  FOP  PLOT  TYPE  D 

1300  CALL  3RF  IXPI*  0*«  XMA*  1  *  *  XNf «  0«*  OX* «05* I PA  II T ) *2* 1 * J 000# T 1 1  ) « 
IFNI 11*  1*  1*  2H75.-1.  0) 

1310  CONTINUE 
1320  CONTINUE 
1330  IRUN  ■  I  PUN  +  1 

IF  ( 1  RUN  -  NRUN)  1340*  1340*  1350 
1340  GO  TO  iLlNE 
1350  CALL  INTVL  12) 

1360  FORNAT  I19H0TIME  IN  SECONDS  ■  .  F10.2) 

PRINT  1360.  Z 
GO  TO  100 
END 

T  SUBTYPE* DATA 
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